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ABSTRACT 

We present a determination of the local (z ~ 0) luminosity function of optically selected type 1 (broad-line) Active Galactic Nuclei. 
Our primary resource is the Hamburg/ESO Survey (HES), which provides a well-defined sample of more than 300 optically bright 
AGN with redshifts z < 0.3 and blue magnitudes B S 17.5. AGN luminosities were estimated in two ways, always taking care to 
minimise photometric biases due to host galaxy light contamination. Firstly, we measured broad-band Bj (blue) magnitudes of the 
objects over small apertures of the size of the seeing disk. Secondly, we extracted Ho- and Hj8 broad emission line luminosities from 
the spectra which should be entirely free of any starlight contribution. Within the luminosity range covered by the HES (-19 i Mb, ~ 
-26), the two measures are tightly correlated. The resulting AGN luminosity function (AGNLF) is consistent with a single power 
law, also when considering the effects of number density evolution within the narrow redshift range. We compared our AGNLF with 
the Ha luminosity function of lower luminosity Seyfert 1 galaxies by Hao et al. (2005) and found a smooth transition between both, 
with excellent agreement in the overlapping region. From the combination of HES and SDSS samples we constructed a single local 
AGNLF spanning more than 4 orders of magnitude in luminosity. It shows only mild curvature which can be well described as a 
double power law with slope indices of -2.0 for the faint end and -2.8 for the bright end. We predicted the local AGNLF in the 
soft X-ray domain and compared this to recent literature data. The quality of the match depends strongly on the adopted translation 
of optical to X-ray luminosities and is best for an approximately constant optical/X-ray ratio. We also compared the local AGNLF 
with results obtained at higher redshifts and find strong evidence for luminosity-dependent evolution, in the sense that AGN with 
luminosities around Mb - -19 are as common in the local universe as they were at z = 1.5. This supports the 'AGN downsizing' 
picture first found from X-ray selected AGN samples. 
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1. Introduction 

A good knowledge of the AGN luminosity function (AGNLF) 
is important for our understanding of the AGN population and 
its evolution, as well as for gaining insight int o the history of 
black hole growth a nd galaxy evolution (e.g. lYu & Tremaind 
l2002t iMarconi et al.ll2004) . Thanks to recent heroic quasar sur- 
veys such as the 2dF QSO Redshift Survey (2QZ) and the 
Sloan Digital Sky Survey (SDSS), large samples of AGN be- 
came available, and the optical AGN/QSO luminosity func- 
tion is today well established over a wide range in redshifts 
(0.3 



5) and luminos i ties flBovle et al 



2004t [Richards et alJl2006t iBongiorno et al 



200C; 



2007 



Croom et al 



Croometal 



2009b^. However, neither of these surveys reaches below red- 
shifts of z ~ 0.3-0.4, because the imposed colour selection cri- 
teria, and also the discrimination against extended sources in 
the 2QZ sample, exclude a large fraction of lower luminosity, 
'Seyfert-type' AGN. The local AGNLF is therefore much less 
constrained than at higher redshifts, despite of its importance as 
a zero-point for studies on the AGNLF and quasar evolution. 

Approaching this problem from the other end, large spec- 
troscopic galaxy surveys are a powerfu l way to unravel the 
AGN content of local galaxy samples (e . g . iHuchra & Burgll992t 
lUlvestad & H"oll200lHHao etai ] 12005 a'). When constructing an 
AGNLF from galaxy surveys, care has to be taken in defining 
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AGN luminosity': Simply taking the optical galaxy magnitude 
will lead to an ill-defined mix of host galaxy and AGN contri- 
butions. It is much better to base the AGNLF on the conspicu- 
ous broad Balmer emission lines which avoids the host galaxy 
contamina t ion pro blem. This approach was re cently adopted by 
iHao et alJ (l2005ah and iGreene & Hoi (|2007|) . But since AGN 
are scarce among galaxies unless the selection i s sensitive to 
very low levels of nuclear activity dHo et al.ll 1997b . any general 
galaxy survey has a low yield, and the samples are dominated by 
weak Seyfert nuclei and rarely reach quasar-like luminosities. 

In this paper we study the population statistics of low- 
redshift quasars and moderately lumin ous Seyferts. We exploit 
the Hamburg/ESO Survey (HES; e.g. I Wisotzki et al . 2000), a 
survey which was specifically designed to address the selection 
of low-redshift AGN with visible host galaxies, and to provide a 
well-defined, complete, and unbiased sample of the local AGN 
population. Here, the term 'complete' is meant in a methodolog- 
ical sense, implying that all AGN selected by the criteria are in- 
cluded in the sample. Since obscured or heavily reddened AGN 
are typically not found in the HES, our sample can only claim 
to be representative of the unobscured 'type 1' AGN with broad 
emission lines in their spectra. We refrain in this paper from re- 
peatedly recalling this fact, but it should be understood that we 
always imply this restriction when using the term 'AGN' . 

This is the third paper investigating the local AGNLF in the 
HES. Koehler et al. ( 1997) used a small sample of 27 objects ob- 
tained during the initial period of the survey to constru ct the first 
published local LF of quasars and Seyfert 1 galaxies. IWisotzkil 
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( l2000bl) exploited ~ 50 % of the HES to study the evolution of 
the quasar luminosity function and provided a re-determination 
of the local LF as a by-product. With the completed HES cov- 
ering almost 7000 deg^ of effective area, the AGN samples are 
doubled in size, as well as augmented by spectroscopic material. 
Thus a new effort is well justified. In the present paper we go 
substantially beyond our previous work not only in quantitative, 
but also in qualitative terms. We first present our data material 
and our treatment of the spectra. After constructing the standard 
broad-band AGNLF we investigate the Ho- and H/3 emission line 
luminosity functions. This enab les us to compare and combine 
our data with the recent work by Hao et al] (2005a) based on the 
SDSS galaxy sample, to obtain a local AGNLF covering several 
orders of magnitude in luminosities. We discuss our results in 
the context of other surveys and luminosity functions. 

Thoughout this paper we assume a concordance cosmol- 
ogy, with a Hubble constant of //q - 70 km s"' Mpc and 
cosmological densi ty parameters Q.„i - 0.3 and Q.^ = 0.7 
dSpergel et al.ll2003l) . 

2. Data 

2.1. The Hamburg/ESO Survey 

The Hamburg/ESO Survey (HES) is a wide-angle survey for 
bright QSOs and other rare objects in the southern hemisphere, 
utilising photographic objective prism plates taken with the ESO 
1 m Schmidt telescope on La Silla. Plates in 380 different fields 
were obtained and subsequently digitised at Hamburg, followed 
by a fully automated extraction of the slitless spectra. In to- 
tal, the HES covers a formal area of -9500 deg^ in the sky. 
For each of the ~ 10^ objects extracted, spectral information 
at 10-15 A resolution in the range 3200 A <, A <, 5200 A is 
recorded. Deta ils about the survey procedure are provided by 
IWisotzki et all ( l2000l) . 

The relatively rich information content in the HES slitless 
spectra enabled us to apply a multitude of selection criteria, de- 
pending on the object type in question. AGN can be easily recog- 
nised from their peculiar spectral energy distributions if they 
contain a sufficiently prominent nonstellar nucleus, i.e. 'type 1' 
AGN. The HES picks up quasars with B $ 17.5 at redshifts of up 
to z ^ 3.2. Several precautions ensured that low-redshift, low- 
luminosity AGN are not systematically missed. For example, 
both the extraction of spectra and the criteria to select AGN can- 
didates contained a special treatment of extended sources, mak- 
ing the selection less sensitive to the masking of AGN by their 
host galaxies. This property makes the HES unique among opti- 
cal survey in that it targets almost the entire local (low-redshift) 
population of type 1 AGN, from the most luminous quasars to 
relatively feeble low-luminosity Seyfert galaxies. It is this prop- 
erty which we exploit in the present paper. 

2.2. Photometry 

Photometric calibration in the HES is a two-step process. We 
first measured internal isophotal broad-band Z?J3 magnitudes 
in the Digitized Sky Survey (DSS) and calibrated these against 
external CCD photometric sequences. At the start of the HES 
around 1990, only the Guide Star Photometric Catalog (GSPC 

' Bj is a roughly rectangular passband defined by a red cutoff at 
5400 A (as imposed by the blue-sensitive photographic emulsion Kodak 
Illa-J) and a blue cutoff at 3950 A. For AGN-like spectra, Bj + 0.1, 
in the Vega system. 




Bj [mag] 

Fig. 1. Effective area of the Hamburg/ESO survey as a function 
of the Galactic extinction corrected B j magnitude. 



iLasker et al.lll988|) was available which typically provided stan- 
dard stars with 8 ^ B $ 14, clearly insufficient to obtain rea- 
sonably accurate photometry near B ~ 17 ... 18, the faint limit 
of the HES. We therefore launched a campaign to obtain our 
own deeper photometric CCD sequences in as many fields as 
possible; the results of this campaign will soon be made pub- 
lic (Schulze et al., in preparation). Augmenting these sequences 
by incorp orating also the Guide Star Photometric Catalog II 
(GSPC-II; 'Buccia relli et"ani2001h . we arrived at a satisfactory 
photometric calibration of all 380 HES fields. 

However, these isophotal DSS magnitudes suffer from a 
number of drawbacks, the most undesirable one (for our pur- 
poses) being the fact that for extended objects, the nuclear AGN 
contribution tends to be drowned out by the host galaxy. Using 
isophotal magnitudes would have a detrimental effect on the 
estimation of a proper AGN luminosity function especially at 
the low-luminosity end. Other negative properties of the DSS 
magnitudes include the effects of relatively poor seeing in many 
fields, boosting the number of undesired blends. 

In order to approximate the concept of nuclear magnitudes, 
we optimised the extraction procedure of slitless spectra in the 
HES digitised data to assume either true point sources or point 
sources embedded in a diffuse envelope. In other words, the 
spectra used for the HES candidate selection always refer to an 
area of the size of the central seeing disk only (note also that 
the HES spectral plates were typically obtained under much bet- 
ter seeing conditions than the DSS). We then measured 'nuclear 
magnitudes' by integrating the spectra over the B j passband, and 
calibrated these magnitudes against the DSS using hundreds of 
stars in each field. We demonstrate below (§|4| that this proce- 
dure indeed produces measurements that in reasonably good ap- 
proximation can be taken as basis for nuclear luminosities (see 
also Figs.|5]and|6]below). 

An additional advantage of our approach lies in the fact that 
we thus entirely avoid variability bias: Selection criteria and 
fluxes in each field are all defined for a single dataset. 

For the purpose of the present paper, a good knowledge of 
the survey selection function, and thus of the accurate flux limits 
is important. In the case of the HES, the flux limit varies consid- 
erably between individual fields, mostly due to changes in the 
seeing and night-sky conditions, but also because of plate qual- 
ity. Thus, each field maintains its own flux limit, always corre- 
sponding to the same limiting signal to noise ratio in the slitless 
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spectra. As long as redshift-dependent selection effects are ne- 
glected, the survey selection function is identical to the 'effective 
area' of the survey, combining the 380 flux limits into single ar- 
ray that provides the total survey area as a function of magnitude. 
The effective area of the entire HES is shown in Fig.[T] Notice 
that this array also fully accounts for the losses due to overlap- 
ping spectra and images (see the full discussion in lWisotzki et al.l 
EOOO). 

All magnitudes used in this paper ha ve been corrected for 
Galactic extinction, using the dust maps ofSchlegeletaDOHi 
(averaged ove r each HES field), and the extinction law of 
ICardelU etall ( 1989). Note that for the HES we had previously 
followed a different extinction correction recipe based on mea- 
sured column dens ities of Galactic ne utral hydrogen converted 
into extinction (see lWisotzki et al]|200 0). The main result of this 
change is a slight decrease of the average adopted extinction 
along most lines of sight. 

Figure[T]shows that the average extinction-corrected limiting 
magnitude of the HES is By < 17.3, but with a dispersion of 
0.5 mag between individual fields. There is essentially no bright 
limit, with 3C 273 recovered as the brightest AGN in the sample. 

2.3. Spectroscopic data 

All AGN and QSO candidates brighter than the flux limit in a 
given field were subjected to follow-up spectroscopy for confir- 
mation and accurate redshift determination, altogether approx- 
imately 2000 targets. These observations were carried out dur- 
ing 23 observing campaigns between 1990 and 2000, using var- 
ious telescopes and instruments at the ESO La Silla observatory. 
Known AGN that were recovered by the HES selection crite- 
ria were initially not part of the follow-up scheme. However, 
such objects were often included as backup targets, so that our 
spectroscopic coverage is almost complete for the entire AGN 
sample. For some quasars in areas overlapping with the Sloan 
Digital Sky Survey we could later add spectra from the public 
SDSS database. 

Because of the diversity of telescopes used, the quality of the 
spectra varied significantly, in signal to noise as well as spectral 
resolution. We estimated the latter quantity from the width of the 
night sky O i emission line at 5577 A. A list of individual observ- 
ing campaigns with their estimated spectral resolution and the 
number of AGN from our sample observed is listed in Table [T] 

Although a formal spectrophotometric calibration was avail- 
able for all campaigns, the combined dataset is much too hetero- 
geneous to allow for any consistent direct measurement of emis- 
sion line fluxes. We therefore adjusted all spectra to a common, 
homogeneous flux scale by computing a synthetic B j magnitude 
from each spectrum and matching this to the 'nuclear magni- 
tude' photometry of the HES. This step again also ensured that 
AGN variability is of no importance. 

2.4. The sample 

For the present investigation of the 'local' AGN population we 
selected all AGN from the final HES catalogue (Wisotzki et al., 
in prep.) that belong to the 'complete sample' and that are lo- 
cated at redshifts z < 0.3. This sample contains 329 type-1 AGN. 
For 5 of them we could not obtain spectra of sufficient quality. 
Thus our sample is 324/329 w 98.5 % complete in terms of 
spectroscopic coverage. 

The HES nuclear fluxes were c onverted into abso lute Mb, 
magnitudes using the K correction of lWisotzkil(l2000ah . Figure|2] 



Table 1. Observing campaigns for follow-up spectroscopy 



Date 


Telescope 


No. of Spectra 


Resolution /1/A/i 






4 


1020 


Apr 1991 


3.6 m 


1 


830 


Feh 1 QQT 


2.2 m 


4 


370 


Feb 199? 


3.6 m 


5 


900 


Cpr, 1QQ9 


3.6 m 


4 


700 


Mar 199^ 


3.6 m 


5 


770 


Feb 1994 


1.5 m 


\ 


1900 


Sen 1994 


1.5 m 


2 


1980 


Sen 1 994 

ijv^iJ 1 y y~ 


3.6 m 


4 


680 


Nov 1994 

i 'Hj V 1 yy^ 


1.5 m 


\\ 


770 


Oct 1995 


1.5 m 


30 


650 


Oct 1995 


2.2 m 


\ 


300 


Mar 1 996 

i.v±cii 1. y y\j 


1.5 m 


19 


590 


WCl Vyyvj 


1 ,J 111 


e 

o 




Feb 1997 


1.5 m 


24 


600 


Oct 1997 


1.5 m 


19 


540 


Dec 1997 


1.5 m 


33 


540 


Sep 1998 


1.5 m 


26 


670 


Nov 1998 


1.5 m 


23 


600 


Sep 1999 


1.5 m 


28 


580 


Mar 2000 


1.5 m 


15 


540 


Sep 2000 


1.5 m 


36 


640 


Nov 2000 


1.5 m 


22 


790 




^ I I I \ \ \ \ \ I I I \ \ \ 

0.1 0.2 0.3 

redshift 



Fig. 2. Redshift against absolute magnitude Mb,. The dashed 
line indicates a constant apparent magnitude of 17.5 mag, ap- 
proximately the flux limit of the HES. 

shows the distribution of the sample over redshifts and absolute 
magnitudes. A wide range of nuclear luminosities is covered, 
ranging from bright quasars with Mb, - -26 to low-luminosity 
Seyfert 1 galaxies of only Mb, - -18. 

3. Emission line properties 

3.7. Fitting procedure 

As our sample is defined by z < 0.3, all objects have and 
most also have Ho- visible in their spectra. We measured lu- 
minosities and line widths of the Ho- and H/3 emission lines 
by fitting the spectral region around each line with a multi- 
component Gaussian model. Over this short wavelength range 
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we approximated the underlying continuum as a straight Hne. 
While the structure of the broad lines is generally complicated, it 
has repeatedly been demo nstrated (e.g. Steidel & Sargent 1991; 
iHo et al.|[T997t iHao et al.i r2005b) that in many cases a double 
Gaussian provides an acceptable fit to the BLR lines. We al- 
ways started with a single Gaussian per line. A second (in a 
few cases also a third) component was added only if a single 
Gaussian yielded a poor fit and the addition of a component 
significantly improved it. At our limited spectral resolution, any 
possibly present narrow component of a broad line such as H/3 
is difficult to detect. We did not include an unresolved narrow 
component by default and added such a component only when it 
was clearly demanded by the data. 

The Fell emission complex affecting the red wing of H/3 
can be sufficiently approximated by a double Gaussian at wave- 
length /1/14924, 5018 A for our data. We fixed their positions 
and also fixed their intensity ratio /15018//14924 to 1.28, typical 
for BLR conditions. Thus only two parameters were allowed to 
vary, the amplitude and the line width. The [O iii] AA4959, 5007 
A lines were also fitted by a double Gaussian with the intensity 
ratio fi xed to the theoretical value of 2.98 (e.g. lDimitriievic et al.l 
I2007h . The relative wavelengths of the doublet were fixed as 
well, but the position of the doublet relative to H/3 was allowed 
to vary as a whole. 

For the Ha line complex, the contributions of [Nii], and 
sometimes also [S ii] needed to be taken into account. The latter 
was mostly well separated and could be ignored, but when re- 
quired we modelled it as a double Gaussian with fixed positions 
and same width for both lines. For fitting the [N ii] /l/i6548, 6583 
A doublet we left only the amphtude f ree. The positi ons were 
fixed, the intensity ratio was set to 2.96 dHo et al.lll997h . and the 
line width was fixed to the width of the narrow [O iii] lines. 

With this set of constraints, each spectrum was fitted with 
a multi-Gaussian plus continuum model. We decided manually 
which model fits best, neglecting contributions by Fe ii, [N ii] and 
[S ii] lines if not clearly present. In 6 objects we detected Ho- but 
the S/N was too low to trace HyS. On the other hand, 21 of our 
spectra did not reach sufficiently into the red to cover Ha; further 
8 spectra were too heavily contaminated by telluric absorption 
to produce a reliable Ha fit. For these objects, HyS was readily 
detected. 

Of the 324 spectra of the sample, we thus could obtain rea- 
sonable fits for 318 objects in Hy6, and for 295 in Ha. Figure [3] 
shows some example results, illustrating the range of signal to 
noise ratios and resolutions of the spectral material. 

3.2. Line Fluxes 

From the fitted model we determined the emission line fluxes of 
Ha and HyS as well as the continuum flux at 5100 A. As said 
above, a narrow component of the Balmer line was only sub- 
tracted if clearly identified in the fit. This happened in 46 cases 
for Ha and in 34 instances for HyS. We also measured the line 
widths of Ha and HyS, which were then used to estimate the 
black hole masses for the sample. These results are presented 
in a companion paper (Schulze & Wisotzki, in prep.). 

In order to estimate realistic errors we constructed artificial 
spectra for each object, using the fitted model and Gaussian ran- 
dom noise corresponding to the measured S/N. We used 500 re- 
alizations for each spectrum. We fitted these artificial spectra as 
described above, fitting the line and the continuum and measured 
the line widths and the line flux. The error of these properties was 
then simply taken as the dispersion between the various realiza- 
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Fig. 3. Examples of fits to the spectra, illustrating the quality of 
the spectroscopic data. The line complex of H/3 is shown on the 
left side, the corresponding complex around Ha is shown on the 
right side. The data are represented by black lines, the multi- 
component Gaussian fits to the Balmer lines are shown in red, 
other lines are shown green, and the combined model is over- 
plotted in blue. Each panel also shows the fit residuals at the 
bottom. 

tions. Note that this method provides only a formal error, tak- 
ing into account fitting uncertainties caused by the noise. Other 
sources of error may include: A residual Fe ii contribution; an 
intrinsic deviation of the line profile from our multi-Gaussian 
model, as well as uncertainties in the setting of the continuum 
level. 

3.3. Relation between H/3 and Ha fluxes 

In order to investigate the distribution of AGN emission line 
luminosities we focus on Ha and Hp as the two most promi- 
nent recombination lines in our spectra. It is of interest to look 
at the relation between these two lines. While recent pubhshed 
work on this subjec t has mostly relied on Ha d Hao et ani2005at 
iGreene & Holl2005h . extending similar studies to nonzero red- 
shifts is easier using Hy6. In Fig. |4] we plot the fluxes of the 
two broad lines against each other. As expected. Ha and HyS are 
strongly correlated. To quantify this, we applied a linear regres- 
sion between Ha and H/3 in logarithmic units, using the FITEXY 
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log (fH/3 [erg s 1 cm 2 ]) 



Fig. 4. CoiTelation between the fluxes /(Hy6) and /(Ha). The 
solid line shows the regression result using the FITEXY method. 
The dashed line shows the best fit line with slope of unity. 

method (iPress et aT l [T992h which allows for errors in both coor- 
dinates. We account fo r intrinsic scatter in the relation following 
iTremaine et alj (|2002|) by increasing the uncertainties until a 
per degree of freedom of unity is obtained. The best-fit relation 
found for the line fluxes is 

logiha) = (1.14 ± 0.02) log (/us) + (0.46 ± 0.01) (1) 

where the fluxes are given in 10"'^ erg s"' cm"^. The rms scatter 
around the best fit is 0.15 dex. The relation between Ha and 
Hj6 line flux using the FITEXY method is shown as the solid 
line in Fig. |4] However, a relation with a slope of unity is also 
consistent with the data (with a scatter of 0. 14 dex). When fixing 
the slope to one, the normalisation corresponds to the relation 
fiia = 2.96 /ha, consistent wit h the Case B recombination value 
of 3.1 (e.g. lOsterbrocgil989l) . We conclude that the measured 
Ha and Hfi broad line fluxes give consistent results. 

4. AGN luminosities 

A long-standing issue for the determination of the local AGN 
luminosity function is the problem of how to disentangle nu- 
clear AGN and host galaxy luminosities. Using total or isopho- 
tal photometric measurements will inevitably lead to luminosity- 
and redshift-dependent biases. For example, the Seyfert lu- 
minosity functions d e termined in some ea rlier studies (e.g. 
iHuchra & BuTa 119921: lUlvestad & Hoi |2001|) clearly reflected 
more the distribution of host galaxy luminosities than AGN 
properties. Ideally, AGN and host light should be properly de- 
composed object by object; but as this would require high an- 
gular resolution data for all objects in a sample, such a route is 
presently not possible. 

As a simplified approach to tackle this problem, we intro- 
duced in Sect. |2.2| our concept of 'nuclear magnitudes' measured 
in the HES spectral plates. We did not subtract any host galaxy 
contribution, but we kept it to a minimum by measuring only the 



flux contribution of a nuclear point source. In Fig. |5]we com- 
pare these 'nuclear' magnitudes with the more standard isopho- 
tal measurements in the DSS direct images. While for high- 
luminosity quasars (and for all quasars at higher redshifts, not 
shown in the figure), these two magnitudes give completely con- 
sistent results, there is an obvious discrepancy which increases 
towards lower luminosities. Clearly, the isophotal magnitudes 
are biased for almost all AGN with nuclear magnitudes fainter 
than -23, and useless for low-luminosity Seyfert galaxies. 

Hence it appears that the HES magnitudes are better esti- 
mates of nuclear luminosities than standard isophotal ones; but 
are they good enough? To investigate this we now compare the 
HES magnitudes with the luminosities of the broad emission 
lines. As Ha and H/3 are pure recombination lines, their l umi- 
nositi es should be proportional to the UV continuum (e.g. lYed 
1980). In an AGN spectrum, the fluxes of the broad lines are 
usually conspicuous and can be reasonably well measured even 
when the host galaxy contribution is strong or dominant. We thus 
adopt the Balmer emission line luminosities as proxies for the 
UV continuum luminosity of the AGN, without any host contri- 
bution. 

In Fig. |6]we compare the HES broad band nuclear B j mag- 
nitudes with the luminosities in both Balmer lines. The correla- 
tion is excellent and, more importantly, it extends over the entire 
range of luminosities. To guide the eye, the dashed lines in Fig. |6] 
represent fixed-slope relations Mgj oc -2.5 log L. We see that the 
data come very close to a slope of unity in the case of Mgj vs. 
Ha, while the relation is slightly shallower for Hfi. 

If the HES magnitudes were systematically affected by host 
galaxy contributions we should expect to see a saturation of Mgj 
values at small L. In fact no clear such trend is visible in the 
data, except maybe a small excess of a few points above the lin- 
ear relation in the lower left corner of the right-hand panel. If 
at all, these few objects appear to be the only ones significantly 
affected by host galaxy light when described by the HES mag- 
nitudes. Overall we conclude that the broad band photometry of 
the HES is, in good approximation, a measure of the pure AGN 
luminosities. 
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Fig. 6. Relation between absolute magnitudes in the Bj band and line luminosities (left: Ha, right: H^S). A fixed-slope relation 
Mbj oc -2.5 log L is shown as a dashed line in each panel, for illustration purposes only. 



5. Luminosity functions 

5.1. Luminosity function parameterisation 

The AGN luminosity function (LF) <p{L) is defined as the num- 
ber of AGN per unit volume, per unit luminosity. The number 
of AGN per unit volume and per unit logarithmic luminosity 
<&(L) = dWlogioi) is given by (1)(L) = (L/ logjo 
where is the cumulative luminosity function. 

In the following we present the results in two different ways. 
We first estimate the LF in discrete luminosity bins, expressing it 
in the logarithmic form '^(L) (or in the equivalent form in mag- 
nitudes). We then show the results of fitting these binned LFs 
with simple parametric expressions. As usual, the fit parameters 
are always expressed in terms of the non-logarithmic form (p{L). 

The most frequently adopted parametric form for the AGN 
luminosity function is a double power law: 



<P(L) 



{Liuy- + {Liuy 



(2) 



where L, is a characteristic break luminosity, (p* the normalisa- 
tion and a and fi are the two slopes. 

It will be seen that the local AGN LF is close to a single 
power law, so we will also consider that even simpler form: 



(3) 
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Fig. 7. The differential broad band quasar luminosity function of 
the HES for 0.01 < z < 0.3. The error bars are based on Poisson 
statistics. The arrow indicates a bin with only a single object. The 
black filled symbols show the luminosity function corrected for 
evolution using a simple PDF model; red open symbols show the 
LF not corrected for evolution. The dashed line shows a single 
power law fit to the data. 



Expressed in absolute magnitudes these functions have fol- 
lowing form: 



a)(M) 



JQ0.4(l+a)(M-M.) ^ JQ0.4(1+^)(M-M,) 

for the double power law, and 

(D(M) = Qi'iQ-^Mi+'^m-M.) 



(4) 



(5) 



for the single power law, with O* = 0.4 0* In(lO) for both func- 
tions. 



5.2. Broad band Luminosity Function 

We first present the broad band luminosity function, using the 
full sample of 329 ty pe-1 AGN with z < 0.3. This updates 
the previous work of ,Wisotzkil (l2000bl) . with several improve- 
ments: (i) The sample size has doubled, (ii) The quality of the 
external photometric calibration is improved, (iii) The correc- 
tion for Galactic extinction has been updated, (iv) We now also 
include the effects of differential evolution within the redshift 
range < z < 0.3 (see below). 
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We determined the bi nned luminosi ty function using the 
classical V/Vmax estimator (' Schmidtll968h . The luminosity func- 
tion is then calculated by: 



M+AM/2 



a)(M) 



AM ^ 



(6) 



where AM is the bin size and V^ax is the survey volume in which 
the object k could have been detected within the flux limit of the 
survey and the given redshift range. Recall that the dispersion 
in limiting magnitudes over the 380 HES fields is taken into ac- 
count by the magnitude dependence of the effective survey area 
f^eff (Fig.[B- Thus Vtnax IS givcn by 



r»» dv 

Vmax = J D.ef[{m)—dz ■ 



(7) 



In adopting this form we explicitly assume that the probability 
of finding an AGN in the HES is independent of redshift. While 
this would be too strong an assumption for the full quasar sam- 
ple, it is certainly justified for the restricted low-redshift range 
z < 0.3 considered here. The SEDs of typical type 1 AGN (with- 
out host galaxy contributions) at such low redshifts are distinctly 
blue in the optic al/UV, and significant marked variation with red- 
shift occur beyond z > 0.5. (Note that we do not consider here 
the role of 'red' quasars which would be lost in the HES al- 
together.) If S(m,z) is the redshift- and magnitude-dependent 
survey selection function, we can safely marginalise over red- 
shifts and set S(m,z) Qeff(m). We reiterate that the effec- 
tive area Qett(m) already fully accounts for the 'target sampling' 
incompleteness due to photometric losses, overlapping spectra 
etc. which h as been carefully dete rmined in the same way as 
described by IWisotzki et al.l (l2000h . Furthermore, as essentially 
100% of the follow-up spectra of the HES could be classified 
with high fidelity, we assume that the probability of selecting an 
AGN within the targeted redshift range and brighter than the flux 
limit is unity. 

The resulting AGN luminosity function is shown in Fig. |2l 
covering the range -26 < Mgj < -18 in bins of 0.5 mag. It 
shows remarkably little structure and rises steadily up to the 
faintest luminosities in the sample. At Mgj S -19 there ap- 
pears to be an abrupt break which almost certainly is an artefact, 
indicating the inevitable onset of severe incompleteness in the 
HES sample for very low luminosities. For such objects, the host 
galaxy contribution even to the HES nuclear extraction scheme 
will be substantial, modifying the slitless spectra in a way that 
they no longer can be discriminated from normal, inactive galax- 
ies. It is remarkable that this effect plays a role only for the low- 
est luminosity bins; there is no gradual turnover that might indi- 
cate incompleteness already at higher luminosities (alternatively, 
invoking incompleteness would imply an even steeper LF which 
would be inconsistent with other results, see below). 

If we ignore the lowest luminosity bins affected by incom- 
pleteness, the Mbj AGN luminosity function is consistently de- 
scribed by a single power law with slope a - -2.4. Fitting in- 
stead a double power law to all bins results in the same bright- 
end slope and a break at Mbj - -18.75. A double power law fit 
to these data is apparently not physically meaningful. 

The observation that the local AGNLF is perfectly described 
by a single power law is in excellent qualitative and quantita- 
tive agre ement with previous results obtained in the course of 
the HES dKoehler et al.lll997l:IWisotzkill2000bh . It is however in- 
consistent with z - extrapolations of the double power law 
AGNLF obtained at higher redshifts. We will discuss this point 
in section [641 



Table 2. Binned differential luminosity function in the B j band, 
corrected for evolution. is the number of objects contributing 
per bin. The luminosity function is expressed as number density 
per Mpc^ per unit magnitude. 
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The evolution of comoving AGN space densities with red- 
shift is sufficiently fast that there is a noticeable effect even 
within the range < z < 0.3. To derive a truly 'local' 
(z - 0) AGNLF we have to take evolution into account. The 
W^max formalism (ISchmidtlll968 ^ provides a simple but ade- 
quate recipe to do so. If our sample were unaffected by evolution, 
we would expect to find (V/ Vmax) - 0.5. Our measured value is 
(W^max) = 0.54 + 0.02, implying some evolution. To correct the 
z < 0.3 AGNLF to z = 0, we approximate the evolution within 
this small redshift interval as pure density evolution (PDF), i.e. 
0(z) oc ( 1 + zY"^ ■ We varied the density evolution parameter ko 
until we reached (V/Vmax) ~ 0.5. To increase the leverage we 
performed the same exercise for a larger redshift range, includ- 
ing quasars from the HES sample up to z = 0.6. We found that 
the evolution within this redshift interval is well described by a 
PDE model with ko - 5, essentially independent of the exact 
value of the outer redshift boundary. Notice that beyond the very 
local universe, the HES samples only the brightest quasars, and 
our correction for the most part concerns these high-luminosity 
bins only. For a single power law, density and luminosity evo- 
lution are indistinguishable. Therefore our evolution correction 
does not critically depend on the choice of the actually adopted 
model. 

We then applied the parameterised density evolution to the 
HES z < 0.3 sample to recompute the evolution-corrected z-Q 
AGNLF. Note that we used the objects at z > 0.3 only to con- 
strain the density evolution index ko and not for the luminosity 
function. 

We still obtain a relation very close to a single power law, 
which however is slightly steeper than the uncorrected one. The 
best fit power law slope is now a - -2.6. This evolution- 
corrected AGNLF is also shown in Fig.|2] It is provided in tabu- 
lated form in Table |2l 
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Table 3. Binned Ho' and li/3 emission line luminosity functions. 
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Fig. 8. Binned emission line AGN luminosity functions, for Ha 
(red filled circles) and H/3 (blue squares). Open symbols indicate 
incompletely filled bins. The dashed lines show the best fit single 
power law to the Ha and H/3 data, respectively. 



5.3. Emission Line Luminosity Function 

The emission line luminosity function (ELF) - the number of 
AGN per unit volume per unit logarithmic emission line lumi- 
nosity - is given by 



<D(L) 



A log L 



(8) 



The selection of AGN in the HES in this redshift range (and 
up to z ~ 2.5) is exclusively based on continuum SED proper- 
ties and independent of emission line properties (Wisotzki et al. 
12000 ). Thus the selection function is the same as for the broad 
band luminosity function, and the Vmax values needed for the 
determination of the emission line luminosity function are also 
the same as for the Mbj luminosity function, except for a cor- 
rection factor containing the spectroscopic incompleteness. As 
discussed in § 12.41 our sample has an overall spectroscopic com- 
pleteness of 324/329 X 98.5 %. Of the 324 objects with spectra, 
6 do not show Hfi and 29 do not show Ha, due to insufficient 
spectroscopic coverage or atmospheric absorption. This incom- 
pleteness does not affect any particular object types preferen- 
tially, and we assumed the losses to be randomly distributed. We 
therefore adopted the effective survey area as before, multiplied 
by a factor of 318/329 for H/3 and by 295/329 for Ha, respec- 
tively. Unless explicitly stating otherwise, we always refer to the 
evolution-corrected Vmax values, thus providing luminosity func- 
tions valid for exactly z = 0. 

The resulting ELFs for Ha and H/3 are shown in Fig. [8l 
binned into luminosity intervals of 0.25 dex. The binned values 
are given in Table [3] 

Figure [8] displays a very similar behaviour of the ELF when 
compared to the broad-band LF of Fig. [T] It rises nearly as a 
straight line, i.e. as a single power law, until a sharp cutoff at low 
luminosities indicates the onset of sample selection incomplete- 
ness. Fitting power law relations to the data (again excluding 
the obviously incomplete lowest luminosity bins) gives slopes 
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of UHa = -2.28 and a^g = -2.26, respectively. Fitting a double 
power law improves the fit quality only marginally. We conclude 
that a description of the ELF as a single power law, for the lumi- 
nosity range covered by our data, seems most appropriate. 

6. Discussion 

6.1. Low luminosity-AGN: Comparison witii SDSS 

The sharp cutoff in the binned luminosity functions at nuclear 
luminosities Mbj ^ -19 or log Lhu ~ 42 clearly signals the 
onset of incompleteness in our sample. This luminosity approx- 
imately marks the limit where AGN cease to be conspicuous 
in the optical and tend to be masked by their host galaxies. 
However, deep spectroscopic surveys of galaxies have shown 
that the AGN p henomenon persists down to very low levels (e.g. 
IHo et al.lll997l) . In those cases, the only traceable indicator of 
nuclear activity in the optical are the emission lines, thus the 
statistics have to be expressed i n terms of an ELF. This was re- 
cently performed bv lHao et alj ( l2005al hereafter H05), who se- 
lected a set of ~ 1000 Seyfert 1 galaxies from the Sloan Digital 
Sky Survey (main galaxy sample) to measure Ha line luminosi- 
ties and construct the ELF. The redshift range covered in their 
sample is < z < 0.15, and the luminosity range is (lO-'^'^- 
10"^^) erg s H05 found their ELF to be in good agreement 
with several parametric descriptions, including single and dou- 
ble power laws and also a Schechter function. However, differ- 
ences between these forms become manifest only at their highest 
luminosities, for L(Ha) ~ lO"*^ erg s"' . Over much of the lumi- 
nosity range, their data suggest a single power law. 

Fortunately the high-luminosity end of H05 overlaps quite 
well with the low -luminosity end of our Ha ELF, so that a com- 
parison is straightforward. This is shown in Fig.|9]where we plot 
the Seyfert 1 Ha luminosity function of H05 (adapted to our 
cosmology) together with our Ha LF. The transition from one 
dataset to the other is remarkably smooth, if the incomplete low- 
est luminosity bins in the HES sample are ignored. Over the lu- 
minosity range in common, both ELFs are fully consistent with 
each other. 
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Fig. 9. Comparison of the Ha luminosit y function derived in this 
work (black points), with the work of iHao et alJ ([2005 a) (blue 
asterisks). The black dotted line gives the best fit single power 
law to our Ha luminosity function. The red dashed dotted line 
gives the best fit power law to the combined data set, whereas 
the red dashed line gives the best fit double power law. 



Even more remarkably, the shape of the combined Ha LF 
- now covering almost 5 orders of magnitude in luminosity - is 
still close to that of a single power law. Already the fit to the HES 
Ha LF alone is almost consistent with the H05 LF. Combining 
both samples, we find a best-fit power law slope of a ~ -2.2. 
This fit is shown as the dashed-dotted line in Fig.|9] 

However, there is evidence for some curvature in the LF, as 
the H05 LF alone is flatter (a = -2.02) than the HES ELF. This 
is manifest in a 'bulge' around 10"*^ erg s^^ where the space 
density of the combined LF is above the fitted single power 
law. A better fit is obtained by a double power law breaking at 
logL* = 41.67, with a faint-end slope of a = -1.8 describing 
the SDSS data and a bright-end slope of [i - -2.6 describing the 
HES. This mildly curved relation, shown by the dashed line in 
Fig.|9] traces the combined data extremely well. 

Some caveats are in place regarding the combination of the 
two datasets. Although H05 determined the narrow Ha compo- 
nents separately, their broad-line ELF does not have this com- 



ponent subtracted, whereas we tried to remove it. As such a re- 
moval was possible in only a small number of cases, the dif- 
ferent treatment does not make much difference. More relevant 
might be the possibility of a systematic variation of the narrow 
Ha contribution to the total Ha flux with luminosity. However, 
the narrow Ha LF published by H05 has also a slope of -1.8, 
which indicates that there should be no major bias introduced. 

Another methodical difference lies in the fact that we cor- 
rected our LF for evolution, whereas H05 did not. If we assume 
that the most luminous objects of H05 lie close to their high- 
redshift limit of z < 0.15 and adopting our simple PDF recipe, 
then the H05 space densities would have to be corrected down- 
ward by a factor of 1.15^ ~ 2, i.e. by 0.3 dex; this correction 
would rapidly decrease towards lower luminosities. The net ef- 
fect would hardly be visible in Fig.|9] 

We note that recentlv. fOreene & Hol(l2007h (hereafter GH07) 
derived the AGN Ha LF based on a combined sample of about 
9000 broad line AGN from the SDSS. Their space densities are 
considerably below ours, and the two LFs are highly inconsis- 
tent (the same discrepancy exists between the GH07 and H05 re- 
sults). This inconsistency has been traced back to an error made 
by GH07 in the determination of their V/Vmax values (J. Greene, 
private communication). Therefore the luminosity function as 
well as the black hole mass function presented in GH07 are in- 
correct. Removing the error alleviates the discrepancy, and the 
corrected Ha luminosity function for the low redshift AGN sam- 
ple from GH07 is consistent with the HES luminosity function 
presented in this work. 

6.2. The combined local AGN luminosity function 

Given the good agreement of the local HES LF, which traces the 
bright end, and the SDSS Seyfert 1 LF sampling the faint end, 
and having the mentioned caveats in mind, it is justified to com- 
bine both datasets. Our aim is to present a single best-knowledge 
local AGNLF in the optical, covering the broadest possible range 
of luminosities. Because of its robustness against dilution due to 
host galaxy light, we used the Ha luminosity functions of HES 
and SDSS. We adopted a bin size of 0.24 dex (0.6 mag) as an 
integer multiple of the binned values published by H05. We then 
recomputed the Ha ELF from the HES data for the same bins 
and merged the two datasets, with weights provided by the in- 
verse statistical variances. 

For easy comparison with other AGN luminosity functions, 
especially at higher redshifts, we converted this combined LF 
into two common reference systems: (i) absolute magnitudes 
Mb in the standard Johnson B band; and (ii) bolometric lu- 
minosity units. These conversions involve translating the Ha 
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Fig. 10. Local AGN luminosity function constructed as a combination of the SDSS broad line AGN LF from' Hao et alJ (12005a') and 
our Ha HES bright AGN LF, both converted to the B band. The black solid line and dotted line show the double and single power 
law fits, respectively. The open blue circles show the directly determined broad band (Mb,) LF for comparison. 



data into broad band or bolometric fluxes. There is a very tight 
correlation between Ha luminosity and absolute blue magni- 
tude in the HES (Fig. |6]l. The translation relation is Mb - 
-2.1 (logL(lia) - 42) - 20.1. This relation is covered by our 
data down to Mg - -19. In order to incorporate also the lower 
luminosity SDSS data we now make the somewhat unguarded 
step of extrapolating the translation of Lua to Mb towards lower 
L. This certainly introduces additional uncertainties, including 
the possibility that some of the lowest luminosity AGN could 
have very different spectral energy distributions (for example, 
the structure of the accretion disk might change drastically). On 
the other hand, there is no reason to expect such a change to oc- 
cur just at the transition luminosity from HES to SDSS, so some 
degree of extrapolation is most probably justified. 

The resulting combined local AGN LF is shown in Fig. [TOl 
ranging from Mg ^ -25 to Mb ^ -15 (with the above caveat). 
Also shown as a dotted line is the best fit single power law, and 
as a solid line the best fit double power law; the fit parameters 
are provided in TableU] The double power law gives a very good 
overall description, but the departure from a single power law is 
not large, albeit statistically significant. There are some minor 
wiggles in the binned LF that are most probably due to under- 
lying unaccounted for systematics; note however that there is no 
trace of the HES-SDSS intersection (rather: transition region) 
around Mg 20. 

For the convenience of the reader we also provide this LF in 
tabulated form, both in terms of the B band and as a bolometric 
LF. For the latter we adopted the lumin osity-dependent bolomet- 
ric corrections of lMarconi et alj(l2004l) . A separate double power 
law fit to the resulting bolometric LF is also provided in TableU] 
Recall that this bolometric luminosity function is valid only for 
broad-line (type 1 ) AGN, without any accounting for obscura- 
tion. 



6.3. Comparison with X-ray selected sampies 

X-ray surveys have made a great impact on our understanding 
of the AGN population, chiefly through their ability to find low- 
luminosity AGN at all redshifts. However, because of the expen- 
sive spectroscopic follow-up, sample statistics are still moderate 
despite considerable efforts. 

The only dedicated effort to estimate an optical local AGNLF 
from a n X-ray selected sample was published by Londish et aT] 
(120001) . Their bright-end slope of /J = -2.1 (without evolu- 
tion correction) is similar to ours, but they obtained a shallow 
faint end slope of a = -1.1, however with considerable error 
bars. Our new results, in particular in combination with SDSS, 
show clearly that such a flat slope is ruled out and that the local 
AGNLF continues to rise towards very faint luminosities. 

A local {z = 0) luminosity function from an AGN sample 
selected in the soft X -ray (0.5-2 keV) band was presented by 
Ha singer et al.l (l2005l) . Their sample is restricted to unabsorbed 
type 1 AGN and is therefore very comparable to ours. To facil- 
itate a compa rison, we again used the L-dependent bolometric 
corrections of lMarconi et al.l (|2004|) to convert our combined op- 
tical AGNLF int o the soft X-ray domain. In Fig.[TT]we compare 
the z = XLF of lHasinger etM (l2005l) with our prediction. The 
X-ray points are represented by the filled blue triangles, which 
are the b inned estimate i n the r edshift shell z = 0.015 - 0.2 of 
Fig. 7 in iHasinger et al.l (1200 5^ corrected to redshift zero. The 
blue dotted line shows their best-fit luminosity-dependent den- 
sity evolution (LDDE) model for z -0. 

For intermediate and high luminosities, the observed XLF 
and the prediction based on the optical AGNLF are in good 
agreement. However, the two LFs disagree strongly at the faint 
end, with the optical LF predicting more than an order of magni- 
tude higher space densities at given X-ray luminosity compared 
to the directly determined XLF. 
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Table 5. Combined binned local AGNLF, based on the SDSS 
broad line galaxy sample (faint end) and the Hamburg/ESO 
Survey (bright end). The bolometric AGNLF has b e en com puted 
using the bolometric corrections bv lMarconi et alj (|2004|) . 
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Is there an explanation for these discrepancies? One possibil- 
ity might be that the adopted conversion from optical to X-ray 
luminosities has been inadequate. In order to explore this option 
we alternatively tried a constant optical/X-ray luminosity ratio; 
the dashed line in Fig. [TT] shows our double-power law relation 
converted with such a relation. Evidently, that predicted LF is a 
very good match to the bright end of the XLF At the faint end 
the two LFs still disagree, but the disagreement is now much 
less. Thus it appears possible that the luminosity dependence of 
the optical/X-ray ratio is much weaker than usually assumed; if 
so, it would certainly help to reconcile the two luminosity func- 
tions. Another possibility is, of course, incompleteness among 
the fainter objects in the X-ray sample. We reiterate however 
that since both optically and X-ray selected samples contain only 
broad-line AGN, any incompleteness due to obscuration should 
be irrelevant in this context. 



6.4. Comparison with tiigher redshifts: Evidence for 'AGN 
downsizing' 

We now perform a direct comparison of our local AGNLF with 
luminosity functions based on surveys that probe mainly the 
higher redshift AGN population. This is interesting not only be- 
cause a comparison with z = provides the longest leverage 
in redshifts, but also because the local luminosity function can 
be traced to very faint luminosities and high space densities. A 
full investigation of the redshift evolution of the AGNLF is out- 
side of the scope of this paper; we limit ourselves to a simple 
comparison between the z = LF determined above and para- 
metric representations of the AGNLF evaluated at z = 1.5, the 
latter chosen as a representative point at moderately high redshift 
where several surveys have been able to leave their marks. 
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Fig. 11. Comp a rison of the soft X-ray LF of type 1 AGN by 
iHasinger et al.l (l2005h with the prediction based on the optical 
local AGNLF. The filled blue triangles show the binned X-ray 
(0.5-2 keV) LF within z = 0.015 - 0.2, corre cted to z = 
The bl ue dotted line shows the LDDE model of IHasinger et al.l 
(l2005h for z = 0. The filled black circles and the solid line 
show the binned data and the double power law fit to the opti- 
cal HES-(-SDSS LF, converted to soft X-rays using a luminosity 
dependent correction. The dashed line shows the same double 
power law fit, but converted to soft X-rays using a constant cor- 
rection. 
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Fig. 12. Shape evolution and 'downsizing' of the AGNLF be- 
tween z = L5 and z = 0. Filled circles and solid line denote 
our combined HES-i-SDSS AGNLF at z = 0. The blue dashed 
and red dotted lines show the best fit mo del to the data of the 
z = 1.5 AGNLF of t he 2SLA0-HSDSS (iCroom et alJ l2009bh 
and the VVDS-hSDSS jBongiorno et alj2007l) . respectively. The 
thin black dashed, dashed-dotted and dotted lines show the z = 
extrapolations of the LEDE model and the LDDE model by 
Cro om et all ([2009b) and the LDDE model bv lBongiorno et all 
(.20071) . respectively. 
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The results of this comparison are displayed in Fig. [12] The 
datapoints and the solid line show our combined local AGNLF. 
The thick dashed line shows the z — 1.5 fit to the optical 
AGNLF, very recently obtained by ICroom et al.l (l20 09b). com- 
bining QSOs from the 2SLA survey dCroom et aljr2009 a) with 



the results from the SDSS jRichards et al.1 l2006h at the bright 



end. They found their data to be in good overall agreement with 
a double power law LF and use a luminosity and density evolu- 
tion model that also included evolution in the bright-end slope 
(LEDE). 

As a second refe rence we considered the AGNLF by 
iBongiorno et al.l (l2007l) (dotted line), derived from a combina- 
tion of the SDSS at the bright end and the fa int type 1 AGN 
sample from the VVDS dGavignaud et alj|2006h . This sample is 
noteworthy in that it is certainly the most complete set of low- 
luminosity AGN at substantial z, as it is purely flux limited and 
not affected by any colour or morphological preselection. As 
seen in Fig. [12] the two z = 1.5 luminosity functions are highly 
consistent with each other. 

We then formally extrapolated these AGNLF models to z = 
0, shown as thin dotted line for the luminosity depend ent den- 
sity evolution (LDDE) model by Bongiorno et al.l (l2007l) . as thin 
dashed line for th e LEDE model and as dashed dotted line for the 
LDDE model bv lCroom et al.l (l2009bl) . AU three extrapolations 
are clearly not in agreement with the local AGNLF, in several as- 
pects: They all ov erpredict the spac e density of bright AGN. The 
L DDE model by iBongiorno et al.1 ( |2007|) and the LEDE model 
bv lCroom et"aLr(l2009bl) also underpredict the spa ce density of 
low-lum inosity AGN, while the LDDE model by ICroom et al.l 
(I2009bl) is roughly consistent with our data at the faint end. 
These mismatches demonstrate that the local AGNLF provides 
indeed additional and independent constraints for the evolution 
of the AGN luminosity function. The above evolution models 
have been derived from higher redshift data which they fit very 
well; evidently, these models may not be extrapolated outside 
the range where they were observationally established. Thus, a 
direct determination of the AGNLF at all z, including z ~ 0, is 
essential. 

Comparing the z = 1.5 iBongiorno et alj (|2007|) or 



ICroom et al] (l2009bl) and the local AGNLF reveals a striking 
change in the shape of the luminosity function: The pronounced 
break visible at higher z is almost absent in the local LF. In terms 
of space densities, this implies that while high-luminosity AGN 
were much more frequent at high redshifts, AGN with nuclear 
luminosities around Mb ~ -19 are as common in the local uni- 
verse as they were at high z. For somewhat fainter AGN this 
relation might even be reversed, although there are too many un- 
certainties to make such a claim. 

Such a 'downsizing' behaviour of the AGN population 
was first dete c ted th rough X-ray surveys (lUeda et alJ I200I 
iHasinger et al.l 20051) and recently also in optical surveys 
dCroom et al.ll2009bl) . Here we strongly confirm the presence of 
this 'downsizing' behaviour in the optical AGNLF which be- 
comes increasingly prominent at low redshifts. We especially see 
clear evidence for evolution in the faint-end slope of the AGNLF. 
Note that this conclusion does not hinge on our extrapolation of 
the Wo-Mb relation, as the crossing of the z = 1 .5 and z = LFs 
occurs at luminosities still covered by the HES. 

7. Conclusions 

We have presented a new determination of the local (z ~ 0) lumi- 
nosity function of broad-line Active Galactic Nuclei. Our sam- 
ple was drawn from the Hamburg/ESO Survey and contains 329 



quasars and Seyfert 1 galaxies with z < 0.3, selected from sur- 
veying almost 7000 deg^ in the southern sky. As a central fea- 
ture, our broad-band magnitudes were measured in the survey 
data with a point-source matching approach, strongly reducing 
the contribution of host galaxy flux to the inferred AGN lumi- 
nosity. Compared to our previous work we have not only sub- 
stantially increased the statistical basis, but also added a number 
of methodical improvements. 

In the construction of the broad band {B j) luminosity func- 
tion, we now included the effects of differential number density 
evolution within our narrow redshift range, S z < 0.3. Since 
the most luminous AGN tend to be located near the outer edge 
of that range, ignoring evolution makes the luminosity function 
appear slightly too shallow. We find that the evolution-corrected 
local luminosity function within -19 ^ Mb, 5 -26 is well- 
described by a single power law of slope a = -2.6, still signif- 
icantly shallower than the z = extrapolation of the AGNLF 
measured at higher redshifts. 

As a second and independent measure of AGN power we 
investigated the distribution of Balmer emission line luminosi- 
ties, in particular the broad Ho' and Hfi lines. These lines can be 
detected and accurately isolated in optical spectra even of low- 
luminosity AGN where the host galaxy is bright compared to 
the nucleus. We found a very tight correlation between Ha lumi- 
nosities and broad band absolute magnitudes Mbj over the entire 
luminosity range of our sample, confirming that host galaxy con- 
tamination to the Mb, magnitudes is unimportant. 

We constructed the broad emission line luminosity functions 
for Ha and Hfi, and found them to agree well with the broad band 
LF. In particular, there is again no trace of significant curvature 
over the covered luminosity range, and a single power law is still 
sufficient to describe the shape of the LF. 

We found excellent consistency between our data and the Ha 
emission line luminosity fun ction of low-lu minosity AGN de- 
termined from the SDSS by dHao et al.ll2005al) . While the two 
datasets are complementary in luminosity coverage, our low- 
luminosity end overlaps very well with their high-luminosity 
end. The SDSS data seamlessly continue the rise of the LF 
towards the domain of low-luminosity Seyferts. The compar- 
ison with SDSS also delineates clearly that below L(Ha) ~ 
lO'*^ erg s"' (or Mb, ~ -19), the HES sample becomes heav- 
ily incomplete; this we suspected already from the shape of the 
HES luminosity function alone. 

We combined the HES and SDSS results into a single z = 
AGN luminosity function covering more than 4 orders of mag- 
nitude in luminosity. This remedies a long-standing shortcom- 
ing of AGN demographics: Despite the heroic survey efforts, 
there was no really well-determined local luminosity function 
that could serve as anchor for a global take on AGN number 
density evolution. The combined local AGNLF is still amazingly 
close to a single power law, but it definitely shows curvature. A 
good description is provided by a double power law with slopes 
a = -2.0 andyS = -2.8. 

Comparing the combined local AGNLF with determinations 
at higher redshifts, we find strong evidence for luminosity- 
dependent evolution, in the sense that weak AGN experi- 
ence a much weaker number density decline, or no decline 
at all, than powerful quasars. This behaviour is well estab- 
lished from X-ray surveys, where a systematic shift of the 
peak in comoving space density with lumi nosity is observed 
(e.g. lUeda et ani2003t IHasinger et allilOOSl) . Known as 'AGN 
downsizing', it is presumably related to the anti-hierarchical 
mass dependence of black hole growth (e.g. Heckman et al] 
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iMerloni & Heinzll2008l: [Shankar et alj|2009l) . The steepening of 
the faint end slope towards low redshift can be understood in this 
scenario by the change of quasar lifetime with peak luminosity, 
and hence black hole mass. The more massive black hole AGN 
die more quickly than lower mass black hole AGN and are there- 
fore not observable in their decaying stage of their light curve, 
whereas lower black hole mass AGN are observable in their less 
lumino us stage an d contribut e significantly to the faint end of the 
LF ( IHopkins et al. 2006; G avignaud et alj|2008l) . Thus the faint 
end of the luminosity function should consist of a mixture of low 
mass black holes accreting at a high rate and higher mass black 
holes with low accretion rates. Investigating this question will be 
the subject of a forthcoming paper 
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